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Abstract 

Mathematical results on some models describing the motion of a tracer particle 
through a Bose-Einstein condensate are described. In the limit of a very dense, 
very weakly interacting Bose gas and for a very large particle mass, the dynamics 
of the coupled system is determined by classical non-linear Hamiltonian equations 
of motion. The particle's motion exhibits deceleration corresponding to friction 
(with memory) caused by the emission of Cerenkov radiation of gapless modes 
into the gas. 

Precise results are stated and outlines of proofs are presented. Some technical 
details are deferred to forthcoming papers. 



1 Background from Physics 

In this introductory section, we recall some results of an analysis presented in a 
companion paper [3J. 

The physical system studied in this paper consists of a heavy, non-relativistic 
tracer particle interacting with the non-relativistic atoms in a Bose gas. Our 
purpose is to introduce some mathematical models describing this system and to 
analyze their properties in the limiting regime where the density of the Bose gas 
becomes very large, the interactions between atoms in the gas become very weak, 
and the mass of the tracer particle is very large. This regime is commonly called 
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" 'mean- field limif. In this limit, the dynamics of the system is determined by clas- 
sical, nonlinear Hamiltonian equations of motion. If the gas exhibits Bose-Einstein 
condensation these equations describe a process of emission of gapless modes into 
the gas by the tracer particle analogous to the Cerenkov radiation emitted by a fast 
charged particle moving through an optically dense medium, [SJEJI3]. The emission 
of Cerenkov radiation causes a deceleration of the particle's motion. Our models 
thus describe a Hamiltonian mechanism of friction with memory Our aim is to 
describe mathematical results on the simplest model exhibiting this mechanism. 

We start by considering a gas of n bosonic atoms confined to some cubical region, 
A, of physical space M 3 . The pure states of this gas are unit rays in the Hilbert 
space 

Tiff = £ 2 (A,d 3 x)®° n (1.1) 

where ® s denotes the symmetric tensor product expressing our assumption that 
the atoms obey Bose-Einstein statistics. The Hamiltonian of the gas is given by 
the operator 

n 

^:=-E^< + A E <K*--*)> (1-2) 

k=l l<k<l<n 

where m is the mass of an atom, denotes the Laplacian on C 2 (A,d 3 x) with, 
e.g., periodic boundary conditions at <9A, A > is a coupling constant, and is a 
bounded two-body potential of positive type and of fast decay at infinity. 
The state space of the tracer particle is given by 

H p := C 2 (A,d 3 x) (1.3) 

and its Hamiltonian is chosen to be 

H P :=~A^ + V(X) 7 (1.4) 

where M is the mass and X G A the position of the particle, and V(X) is the 
potential of an external force. The interaction of the particle with the gas is given 
by 

n 

Hj-.^gJ^Wixk-X), (1.5) 

k=l 

where W is a bounded two-body potential of fast decay at infinity, and g is a 
coupling constant. 

The Hilbert space of the total system is 

=U p ®H { g\ (1.6) 
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and the total Hamiltonian is the operator 

H<ri = H p + Hg> + Hr, (1.7) 

where H p stands for H p ® and -f^ - " 1 for 7|^ p <g) -ffg^. 

We will be interested in studying the properties of this system in the thermo- 
dynamic limit 

A / R 3 , with p := — kept fixed, (1.8) 

where |A| denotes the volume of A and p is the particle density of the gas. For 
this purpose, it is convenient to describe our model in the formalism of second 
quantization. 

The bosonic Fock space of the Bose gas is the Hilbert space 

To := ©£*?4 B \ (1-9) 
the total Hilbert space is given by 

U = U V ®7 G , (1.10) 

and the Hamiltonian by 

H = ®Z H n \ n(n) . (1.11) 

In order to tune the density of the Bose gas in states of small total energy to a 
desired value ~ p, we replace the operator H by 

H := ®™ =0 (H (n) - pn)\ H n + constant. (1.12) 

Note that, on states of a fixed number, n, of atoms, H differs from H just by 
the constant - pn + constant, and hence describes the same physics. Setting 
p := A0(O)p - |0(O), constant := -gW(0)p + |0(O)p 2 |A|, where A indicates 
Fourier transformation, we find for H, in the formalism of second quantization, 
the expression 

H = -^A x + V(X) + f^(Va*)(x)(Va)(x)d 3 x 

+g fW(x- X)(a*(x)a(x) - p) d 3 x (1.13) 

+ f / J{a*{x)a{x) - p)<f>(x - y){a*{y)a{y) - p) d 3 xd 3 y, 

where reference to the integration domain A is omitted and where the operators 
a*(x) and a(x) are the usual creation- and annihilation operators on Tg- They are 
operator- valued distributions on Tg satisfying the canonical commutation relations 

[a#(x), a*(y)\ = 0, [a(x), a*(y)\ = S(x - y)I, (1.14) 
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a* := a or a*. 

We are interested in analyzing the dynamics of the system described by Equa- 
tions (1.10), (1.13) and ( 1.14) in the thermodynamic limit, A /* IR 3 . In particular, 
we would like to determine the states of lowest total energy, i.e., the ground states 
of H. This poses a formidable mathematical problem that is far from being under- 
stood rigorously. In order to simplify matters, we propose to study the following 
limiting regime ("mean-field limit"): Let iV > be a parameter that will eventually 
be taken to oo. We set 



p = N^, A = N- % X g 2 , M = NM , V(X) = NV {X) 



r 



(1.15) 



where po, g, Xo, M and Vq(X) are fixed (i.e., iV-independent). For physical moti- 
vation underlying the choice (1.15), see [3]. We set 

u#f~\ ._ N~^a*(x) - 



b\[x) 




Po 

2 ' 
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The operators b* N (x) and &/v(x) are creation- and annihilation operators satisfying 
the commutation relations 



[b*(x),b*(y)]=0, [b N (x),b* N (y)] = N- 1 5(x-y)I. (1.17) 
In the new variables the Hamiltonian H introduced in Ql.l3| ) is given by 



H = NH 



N 



(1.18) 



where 
H N 



27V 2 Mi 



+ V (X) + J ±(Vb* N )(x)(Vb N )(x) d 3 



X 



;i.i9) 



+g J W(x - X){b* N (x)b N (x) + yf^(b* N (x) + b* N (x))} d 3 x 

+ ¥ / JK(x)bN(x) + ^(b* N (x) + b* N (x)Mx - y) x 
[b* N (y)b N (y) + ^(b* N (y) + b* N (y))} d 3 xd*y. 

The Schrodinger equation for the time evolution of states, ty t , i n the Hilbert space 
~H introduced in (1.10) has the form 



9 T 



at 



(1.20) 



Considering Equations (1.15), (1.17), (1.19) and (1.20), we see that AT 1 plays 
the role of Planck's constant h and that the mean-field limit, A" — > oo, is equivalent 
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to a classical limit, (h \ 0). From Equations (1.17) through (1.20) it is easy to 
heuristically derive the classical Hamiltonian dynamics emerging in the mean-field 
limit. 

The phase space of the limiting classical system is given by 

T :=M 6 x "H 1 ^ 3 ), (1.21) 

where , H 1 (1R 3 ) is complex Sobolev space over IR 3 . 

We choose the usual position- and momentum coordinates (X, P) on IR 6 and 
complex coordinate functions ((3(x), j3(x)) on 'H 1 (]R 3 ). The standard symplectic 
structure on T yields the Poisson brackets 

{X\ X>} = {P t , P 3 } = 0, {X\ P 3 } = Fj, (1.22) 

= 1, 2, 3, and 

{(3(x), /%)} = {0(x), 0(y)} = 0, {(3(x), 0(y)} = id(x - y). (1.23) 

Corresponding to the Hamiltonian Hn we consider the Hamilton functional 

H = H(X,P;pJ) 

■= £- Q +Vo(X) + J^WP(x)\ 2 d 3 x 

+g j W(x - X){\f3(x)\ 2 + 2^Re(3(x)} d 3 x 

+ ¥ / JM X )\ 2 + 2^ReP(x)]<P(x - y)My)\ 2 + 2^Re/3(y)) d 3 xd*y. 

(1.24) 

Setting a(x) = (3(x) + J~^i the Hamilton functional % is seen to be given by 

« = 2& + ^o(X) + ^/|V«(x)| 2 rf 3 x 

+ gJW(x-X)(\a(x)\ 2 - g) d 3 x (1.25) 

+ Y//(l«(^l 2 -?M*-y)(|a(y)l a -?) rf3 ^ 

and the equations of motion are found to be 

*t = -0, P t = -VV(X t )+gj(VW)(x-X t )(\a t (x)\ 2 -p (1.26) 

and 

ia t (x) = (~A + gW(x - X t ))a t (x) + X g 2 (<p * (\a t (x)\ 2 - ^))(x)a t (x), (1.27) 
2m g £ 

where * indicates a convolution. Setting W = 0, we find the explicit ground state 
solutions 

P t = 0, X t = X* : a minimum of Vq(X), 

(i-28) 
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where 9 is an arbitrary phase. These solutions show that, for p > 0, the continuous 
gauge symmetry 

a(x) — > e a(x), a(x) — > e~ l6 a(x) 



of the Hamilton functional (1.25) is spontaneously broken in the ground states, 



which corresponds to Bose-Einstein condensation. One then expects that the Bose 
gas exhibits gapless (Goldstone) modes at zero temperature. In [3], the ground 
state solutions are constructed for W ^ 0. 

It has been shown in [I] that, in a finite periodic box (torus) A, the dynamics 
of the quantum system with initial conditions given by a "coherent state" in "H 
is given by evolving the coherent state along a solution of (1.26), (1.27), 
harmonic quantum fluctuations of amplitude 0(- 



with 



around the classical solution, 
as N — > oo. The spectrum of these fluctuations is found by linearizing Equations 
(1.26), (1.27) around a given classical solution. For W = 0, and choosing the 
classical solutions to be given by a ground state (1.28), the frequency spectrum 



of the harmonic quantum fluctuations of the Bose gas can be found by passing to 
the "Bogoliubov limit" g — > 0. One then finds the dispersion law 



u(k) = \k\ 



k 2 
4m 2 



rn 



as \k\ —7- 0, 
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where k e A* is the wave vector of a normal mode and t>* = y ~^<^(0) is the speed 
of sound in the Bose gas; see [H E], and [7] for a somewhat different approach. 

It is a challenging open problem to find out whether the mean-field limit — > oo 
and the thermodynamic limit A /* M. 3 can be interchanged. (It is not surprising 
that this problem is difficult, because there is no mathematically rigorous under- 
standing of Bose-Einstein condensation in an interacting Bose gas, in the thermo- 
dynamic limit.) However, for an ideal Bose gas (A = 0) and in the Bogoliubov 
limit (g = 0), this problem is understood, [1]. 

From now on, we suppose that we first pass to the mean-field limit and then to 
the thermodynamic limit or that we consider models in the Bogoliubov limit. 

Our analysis shows that, for vanishing potentials W and Vo, the energy- momentum 
spectrum of the linear quantum fluctuations of the system is as indicated in Figure 

m 



G 



*- IPI 

We observe that, for \P\ > P, = Mf», the energy, of a tracer particle with 
momentum P is embedded in the continuous energy spectrum. Standard ideas of 
resonance theory suggest that if the interactions between the tracer particle and the 
atoms in the Bose gas are turned on, i.e., for W ^ 0, a state of the tracer particle 
corresponding to steady motion at a velocity v = jj, with \v\ > v„, is unstable and 
decays into states of smaller velocity by emission of gapless modes into the Bose 
gas (Cerenkov radiation); i.e., the particle motion undergoes deceleration until the 
speed of the particle has dropped to a value below v*, see [5J |3|. Spectral results 
that go in the direction of confirming this picture have been established in j2]- To 
summarize, we note that, apparently, the motion of a tracer particle interacting 
with the atoms of a condensed Bose gas exhibits friction by emission of Cerenkov 
radiation until the speed of the particle is below the speed of sound of the Bose 
gas. 

In this paper, we present mathematical results confirming this picture for the 
simplest models corresponding to Ao = 0, i.e., for an ideal Bose gas, in the mean- 
field limit. In these models, the speed of sound vanishes, i.e. v* = 0, and the 
friction mechanism is at work until the particle comes to rest. In the variables 
(X, P, (3, (3), the equations of motion are given by 

Y - p t 
A * - Aft' 

P t = -WoPO +gf(VW x *)(x)(\(3 t (x)\ 2 + 2^Ref3 t (x)) d 3 x 

and 

ifr(x) = (h Xt )(3 t (x) + ^- W Xt (x), (1.30) 
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where W x * 



x 



W{x - X) and 



h 



x 



1 

2m 



A + gW 



x 



In order to simplify our notations, we set m — 1 and v /po =: K - We assume that 
Vq is smooth and W is smooth, spherically symmetric, of rapid decay at infinity, 
and such that J W(x) d 3 x = W(0) ^ 0. Moreover, for reasons of stability of the 
system, we choose g to be so small that h x does not have any bound states and 
zero-energy resonances. 

It is reasonable to expect that 



P t -> 0, X t ->• Xoo, 
where Xoo is a (local) minimum of Vo, and 

Pt -> 0, i.e., A -> -/t(/i x -) 



1 W Xo 



1.32) 
1.33) 



as time t tends to oo. We are interested in establishing (1.32) and (1.33) and to 
find out how \Pt\ tends to 0, as t — > 0, for a physically reasonable class of initial 
conditions and sufficiently small values of g. 

In order to develop some intuition into the behavior of \P t \ for large times t, we 
consider a steady motion of the tracer particle under the influence of a constant 
external force, i.e. for Vq(x) = —F ■ X, in the limit where g — > 0; (the "B-model" 
in the classification of [3]). Such a motion is a solution of the following equations: 

X t = v, = P t = F + 2K J {VW Xo+vt ){x)ReP t {x) d 3 x, (1.34) 

where (Xq, Pq = Mqv) are the initial conditions for the tracer particle, and 

f3 t {x) = 7 (a - Xq- vt), 

where 7 solves the equation 

1 



iv ■ (V7)(x) 



: (A7)(x) + kW(x). 



The Fourier transformation, j(k), of the solution of (1.35) is given by 

W(k) 



-2k- 



1.35) 



1.36) 



(k — v) 2 — v 2 — iO 

Note that 7 is singular on the sphere {k\\k — v\ — \v\}. Changing variables, p : 
k — v , and plugging (1.36) into (1.34), we find that 



-Re4Kji(p + v) ] -^±^d s p 

Ik-k J(p + v)\W(p + v)\ 2 5(p 2 -v 2 ) d 3 p. 



1.37) 
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(1.37) can be viewed as an instance of Fermi's Golden Rule; see [3]. If W is of 
short range and \F\ is small, so that \v\ will turn out to be small, too, we may 
replace \W(p + v)\ 2 by 1^(0) | 2 on the right hand side of (1.37). Then 

F « (Atc) 2 kW(0)v 2 v, (1.38) 

where v is the unit vector in the direction of v. Not surprisingly, it follows that F 
and v are parallel. 

Let us now assume that F is taken to 0, but that, at large times when \P t \ is 
very small, with || rs v = constant, the state of the Bose gas is close to (3 t (x) = 
7(x — X — vt), where 7 is the solution of (1.35) just constructed. Then the motion 
of the tracer particle can be found by solving the equations of motion 

M v = -(Atv) 2 kW(0)v 2 v. (1.39) 

Choosing initial conditions P = M u e z (where e z is the unit vector in the z- 
direction), we find that P t = M u t e z , where 

u t = -Cu 2 , (C = (47t) 2 /€M - 1 #(0)) 

so that 



u t = (Ct 



M 



Thus 

and it follows that 



|P t | = O(-), as t -»■ 00, (1.40) 

\X t \= [ u t dt 00, (1.41) 
Jo 

as t — > 00. 

This heuristic analysis also shows that if = 0(\k\~ e ), e > 0, or for a 

particle in a confining external potential Vq with non-degenerate minima, or for a 
two-dimensional system with W(0) 7^ 0, but Vo = 0, X t is predicted to converge 
to a finite point X^, as t — > 00. 

The problem with these arguments is that the initial condition of the Bose gas, 
Po(x) = j(x — Xq), with 7 as in (1.35), (1.36), is highly singular and has infinite 
energy, and this is unnatural, physically. We propose to analyze the motion of the 
particle with initial conditions given by 

(X t=0 , P t =o) = {X ,M v), \v\ small, (3 t=0 = A) 

where /3q is small and of reasonably fast decay at infinity. Then the effective 
equations of motion for the particle exhibit memory effects, and it becomes more 
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subtle to find out how \P t \ decays to 0, as t — > oo. This is the problem we wish 
to solve in this paper for the simplest model with Ao = and g — > 0. We hope to 
treat models with Ao > in the "Bogoliubov limit" , g — > 0, in future work. 



2 Mathematical Results on Friction in an Ideal Bose Gas, 
and Strategy of Proofs 

We set the coupling constant Ao to and assume that the potential W 
(Al) is smooth; 

(A2) decays to exponentially at oo; 

(A3) is spherically symmetric; and 

(A4) has the property that W(0) ^ 0. 

We state our results for the so-called B-model, A = 0, g — > 0, V = 0. The 
equations of motion then have the following form (see (1.30), (1.31)): 

(2.1) 



X t = P t = 2Kj(VW Xt ) Refit 



ifit = -\/\fi t + KW x \ 

with K = ^po", W x (x) = W(x - X). 

Our main result is the following theorem. 

Theorem 2.1. For an arbitrary 5 G (0,5*), where 5* > is some constant, there 
exists an e = e(5) > such that if 

||(a;) 3 /3o||2<e, \P Q \ < e, (2.2) 



with (x) = a/1 + \x\ 2 , then 



\P t \ < const t 2 5 ; as t ->■ oo, (2.3) 



and 

-UxrX: 



lim \\fit- 2K(A)- 1 W Xt \\ 00 = 0. (2.4) 

t— >oo 



Remarks. 



(1) We do not have a simple explicit formula for the exponent 5, and it is not clear 
whether it is universal (independent of initial conditions). 

(2) Assumption (A4), which plays a critical role in our analysis, implies that 
A^ 1 W X ^ £ 2 (IR 3 , d 3 x). In contrast, if Xq is positive we expect that fi t remains 
uniformly square-integrable, as t — )■ oo. 
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Next, we outline the strategy of the proofs of our main results: The equations 
of motion (2.1) are a system of semilinear integro-differential equations. We solve 
the equation for /?., given the particle trajectory X, and plug the result into the 
equation for P.. We derive, in Subsection ^1 below, an effective integro-differential 
equation for P. of the form 



P t = L 1 (P.)(t) + L 2 (P)(t) + N(P.)(t), 



(2.5) 



where L\ and L 2 are linear operators from the space T[ 0j t] := {P s G 1R 3 |0 < s < t} 
of momentum trajectories to M 3 , and N is a non-linear operator from Tp ( t] to M. 3 . 
The operator L\ is of convolution type, i.e., 



HP.)(t) 



ds f(t - s)P s 



where 



f(s) = ZRe{W,e'^W), 



(2.6) 



(2.7) 



P s ds. 



with (•, •) the scalar product on £ 2 ( 
operator L 2 is given by 

L 2 (P.)(t) = f(t) 
Let K t £ R, t £ [0, oo), be the solution of the linear equation 
K t = L 1 (K.)(t) = -f*dsf(t-s)K Sl 

K n = 1. 



d?x) and Z some positive constant. The 



(2.8) 



(2.9) 



By Duhamel's principle, Eq. (2.5) can then be converted into the integral equation 

P t = K t P + [ ds K t _ s L 2 (P.)(s) + [ ds K t _ a N(P.)(s). (2.10) 
JO Jo 

It is easy to derive from the properties of the function / that K t = 0(t~^), as 
t — > oo (see (3.11) below). From this we infer that the terms on the right hand 
side of (2.10) do not decay faster than t~%, as t — > oo. This is disappointing, 
and we propose to show that the leading terms on the right hand side of (2.10) 
cancel each other to yield a decay of 0{t~~>~ 5 ), as t — > oo, for some 5 > 0. This 
cancellation is exhibited as follows: We integrate (2.5) over time to find that 

P t = P + [ dsL 1 (P.)(s)+ [ dsL 2 (P.)(s)+ [ dsN(P.)(s). 
Jo Jo Jo 
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Multiplying this equation by K t and subtracting the resulting equation from (2.10 ), 
we obtain 

-t , T \ , ft 



P t {l-K t ) = -K t £ds L 1 (P.)(s) + fids (K t - a -K t )L 2 (P.)(s) 
+ fids (K t -,-K t )N(P.)(s). 



(2.11) 



The first term on the right hand side of (2.11) does not have an improved decay 
in time, yet. We rewrite the second term so as to cancel the leading contribution 
to the first term. Using ( |2.8 ), we find 

- fids (K t _ s - K t )f(s) fidu P u 

= fi ds (K t ~ s - K t )f(s) fi du P u - fi ds {K t - a ~ K t )f(s) fi du P u ^ ^ 

= fi ds (K t . s - K t )f(s) fi du P u - fi ds K t . s f(s) fi du P u 
+K t fi ds f(s) fi du P u . 

Using expression ( |2.7 ) for f(s), the last term on the right hand side of (2.12) is 
found to be given by 

rt 

2K t Re(W,(iA)- 1 e i ^W) / du P u , (2.13) 

Jo 

where we have used that Re(W, (iA)~ l W) = 0. The first term on the right hand 
side of (2.11) is rewritten as follows 

-K t fids Lx(P.)(s) = -K t fids fiduf{s-u)P u 



-K t fi ds fi du Re (W, e* ^V 1 W)P U 



(2.14) 



= -2K t Re{W, (tAy 1 fi du e i ^r i W)P w 

where we have used integration by parts and Re(W, (zA) _1 iy) = 0. Combining 
( pH2| and gjg we find that 

= 2ZK t Re(W, (iA)- 1 fi ds [e iM ^ - e i ^]W)P s 
+Z fi ds (K t _ s - K t )Re(W 7 e^W) fi du P u 
-Z fids Kt^ReiW^^W) fidu P u 
+ fids (K t _ s -K t )N(P)(s). 



Pt(l - K t ) 



(2.15) 
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Among the terms on the right hand side of (2.15), the third term looks trouble- 



some. In order to show that it has the desired decay, 0(t 
recall that 



with 5 > 0, we 



Z fids K t . s Re(W,e^W) = fids K t _J{s) = L a (#.)(*) = K t 



(2.16) 



Since K t = 0(t 2), we may expect that K t = 0(t 2), and this is indeed the case. 
Making the self-consistent assumption that Pt = 0(t~2~ s ), for some 5 > 0, we 



find that all the terms on the right hand side of (2.15) have the appropriate decay 



in t. In order to extract the decay of P t in t from (2.15) in a rigorous manner 



we would like to convert (2.15) into an integral equation for P t by dividing both 
sides by 1 — K t . Since K t — > 1, as t — > 0, some care is needed for small values 
of t. Since we know that Kt — > 0, as t — > 00, it suffices to "wait long enough" 
before dividing by 1 — K t . This amounts to dividing the time axis [0, 00) into 
two subintervals, [0,T] and [T, 00), where T is chosen so large that K t < e, for 
t>T, with e small enough. On the interval [0,T], we use standard existence and 
uniqueness theorems to solve Eq. (2.10), and, on the interval [T, 00), Eq. (2.11) is 
rewritten appropriately with an inhomogenous part depending on {P s }o< s <t- The 
resulting equation is interpreted as an equation on a Banach space 

B s ,t ■= {P\ \t^ +s \\P t \ uniformly bounded for t G [T, 00)} 



with the norm 

ll-P-lkr := ||X[T,oo)^ 2+5 -P||oo 
where xi denotes the characteristic function of the interval J C R. The idea is then 
to show that, for an appropriate choice of 5 > 0, the equation can be solved in B$,t 
by a standard fixed-point theorem. 

We will discuss some key elements in applying the fixed-point theorem in Sub- 
section 13.21 



3 Some Technical Details of the Proof 



In this section we explain two elements of our proof: one is to derive Eq. (2.5), the 
second one is to show how to make a fixed-point theorem applicable. For a more 
detailed presentation we refer to a forthcoming paper. 
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3.1 Derivation of (|2~5j) 



We begin by recasting equations (2.1) in a more convenient form. We define a 
function 5 t by 

A =: -2/t(-A)- 1 ^ + S t . 
The equation for (3 t implies one for 5 t , and using Duhamel's principle, we find that 

S t = e^/3 + 2Ke^(-A)- 1 ^' - £ f ds ^{-A^P, ■ V X W X ° . 

M Jo 

Plugging this equation into the equation for P t we obtain, after some manipulation, 

P t = 2KRe(V x W x \e 1 ff3 }+AK 2 Re(V x W J e 1 f(-A)- 1 W x °- Xt } 



Mo -He(V,W,/> e'^-A)- 1 ^ • V X W X °~ X <) . (3 ' 1} 

Remark: Here we have used the spherical symmetry of W to cancel terms of the 
form (V^W^, S), where S is spherically symmetric. The fact that such terms vanish 
will be used repeatedly. 

Next, we isolate the terms linear in P. 

Re(V a; iy,e 1 ^(-A)- 1 Vy Xo - Xt ) = Re (V X W ,e'^ (-Ay 1 (W x °~ Xt - W)) 

= ^-Re(V x W,e^(-A)- 1 [ dsP s -V x W) 
M o Jo 

1 /"* 
+ — Re(V a; iy,e 1 ^(-A)- 1 / ds P s ■ V X [W X °- X - - W]) . 
M o Jo 

Proceeding with the other term in (3.1 ) in a similar way, we arrive at the equation 

Pt = 



-^Re {V X W , p ds e'^ 1 (-A)- 1 P s ■ V X W) 
+^Re(V x W,&¥ (-A)- 1 Jo ds P s ■ V X W) 
+N(P)(t), 



(3.2) 



where 

mm) 



2nRe{V x W Xt ,e^l3o) 

+^Re (V x W,e^ (-A)- 1 J? ds P s ■ V X [W X °~ 



x s 



W}) 



w Re {V x W,f Q ds e^{-A)^P s ■ V X [W X ^ - W]) 
2kRe(V x W Xt ,e*^ ) 

+|§Re (V X W , e^ (- A)" 1 J* ds P s ■ V x f° ds, P Sl ■ V x W x °~ x n ) 

+|^Re {V X W , $1 ds e^(-A)- l P s ■ V, £ ds, P S1 ■ V X W X ^). 

(3.3) 
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Using that W is spherically symmetric, we have, for k = 1, 2, 3, 



(d Xk W ,e 1 ^ L (-A)~ 1 P s ■ V X W) = -PW{W,e^W) 



(3.4) 



where P^ k > is the fcth component of P, and (3.2) is thus found to have the form 



(2.5). 



3.2 A Fixed Point Theorem 



To render the analysis after Eq. (2.15) rigorous, we reformulate (2.15) in the 
following form: 

P t = j^L(P.)(t) + J ^J*ds(K t _ s -K t )N(P.)(s) 

= ^-KxiT^P-m + Jo ds (K t _ s - K t )[N(P.)(s) - N( X[0 , T) P.)(s)} 
+ j^-L( X[ o,T)P.)(t) + ^ fids (K t _ s - K t )N( X[ o,T)P.)(s) 
= A( X [T,oo)P.)(t) + nX[o,T)P.)(t) 

(3.5) 

where L is the linear part on the right hand side of (2.15), and the terms on the 



last line are defined to be the first and the second lines after the second equality 
sign. 

To meet the criteria of applicability of a standard fixed-point theorem, we need 
to establish the following results: (1) The term F(x[o,t)P) is small in the Banach 
space Bs,t, for some 5 > and T > 0. This will be shown in (3.7), below. (2) 
For some 5 £ (0, |), the map A(x[T,oo) m ) '■ &8,t &s,t is a contraction if restricted 
to a sufficiently small neighborhood of 0. This is verified in (3.9), below. These 
two facts, together with the observation that A(x[t,oo)0) = 0, obviously make a 
standard fixed-point theorem applicable and hence imply global existence of a 
solution in the space Bs,t- 

In what follows we sketch proofs of (1) and (2). 

In order to see that F(x\o,t)p.) is small, for a suitable choice of initial conditions, 
we observe that it is defined in terms of X{o,T)P- an d the initial condition (3 ; 
moreover if P = and /3q = then P. = 0, which implies that F = 0. From 
a simple proof of local existence of solutions we infer the following proposition. 

Proposition 3.1. For any T > 0, there exists some e = e(T) > such that if 
\Pq\ < e then 



(x}- 3 Poh, 



| P t | < T 2 , for any t <T. 



(3.6) 
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Moreover, for any 8 <%, there is an e(T) > satisfying lim e(T) = such that 

T— >oc 

\\F( X[ o,T)P.)h s , T <e(T)- (3-7) 

Now we turn to our sketch of the proof of (2). Define a function Q : [0, |) — > M. + 

by 



1 ^ 



n(5):=- dr _(i_ r )-2[ (r-2 - r - 6 )+r*- d }. (3i 

7T ' 



1 + (1 - r)i V ^ '' L l - 2<5 V ' 

It is easy to find, by direct computation, that the set {5\5 > 0, Q(8) < 1} is not 
empty. One of our key results is the following theorem. 

Theorem 3.2. If we restrict the domain of definition of A(x\t,oo)') t° a ball in 
the space B$ y x around P. = of sufficiently small radius eg, for some 5 > 0, then 
A(x[t,oo)-) is a contractive map. Specifically if \\P.\\b S t; \\Q-\\b S t — e o then 

\\A(X[T,oo)P.) ~ A(X[t,oo)Q.)\\b s ,t < MS) + e]||P. - Q.\\ BttT (3.9) 
where e = e(|P |, || (^c) 3 /3o 1 1 2 , T) has the property that e(0, 0, 00) = 0. 

In what follows we show how to derive (3.9). We focus our attention on two 
terms on the right hand side of ( |2.15[ ), namely 



ri(x[r,oo)<7.)(*) -=- Z I ds [K{t-s)-K{s)]BB(W,e^W) J da x q slX[ T. 



00) j 



and 



r 2 (xFr,oo)?0(*) : = 2ZK t MW, (tA)- 1 [ ds [e 1 ^ 1 - e^]W} X [T,oo)q s , 

Jo 

where we write q for an arbitrary component of the vector P. By direct compu- 
tation, and using observation (2.16), one can show that all the other terms in A 
can be bounded by a small constant e, (see (3.9)). 
We define two functions Qi, Q 2 '■ (0, \) — > M + by 

ni(*) := n * fdr 1 (l- r )-|[ r -|- r -^], 

{l-2d)n J 1 + (1 — r)a 



and 



1 r 1 , 1 



7T 



fi 2 (<5) := — / dr r (l — r) 2 r 2~ 



l + (l-r 



Obviously 

Q{5) =Q 1 {5) + Q 2 (5). 
A key estimate is given in the following lemma. 
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Lemma 3.3. There exists some e(T) > ; with lim e(T) = 0, such that 



T->oo 



||r fc (x[T,oo)9.)lls 4 .T < Pk + e{T)]\\q\\ BStT , k = 1,2. (3.10) 

Proof. We focus on estimating Fx, (the estimate on IT^ can be derived in an almost 
identical way and hence is omitted). 

Using Eq. (2.9) and Fourier transformation we find, after some manipulations, 
that there exists a constant Ck G K such that 

ZK(t) = — ^ — vr-iH + c K r x + o(rh. (3.11) 

V ^ 4|W(0)| 2 K ' 

By standard techniques 

Re(W,e^W) = -2tt§H \W(0)\ 2 + 0(H ) . 



We define i^, M, Ti to approximate these functions, 

ZK(t) := — J 1 7r"tt-i 

4|W/(0)| 2 

3 3 



M := -2|lU(0)| 2 7r2r2 
f 1 := -Z J ds \K{t -s)- K(s)]M(s) J ds x q Sl [1 - X t(si)} . 



We then find that 



t rt 

3 



|fi| < ^ J ds[(t-s) »-t a]s 2 y tfei |g si | 

< n L / ds [(t-s)-*-r?)s-%(t$- s -s*- s )\\q.\\s, T 
(1 - 2d)n J 

1 /"* i / \_i -I 1 -1/ I_A i-A\ n n 

Changing variables, s = rt, we obtain 

\fi\<t-*- s n 1 (S)\\q.\\ s ,T. (3.12) 

In the following we estimate the remainder \Ti — Til- Observe that K and M 
are good approximations of K and Re(W, e % ^W) only when t is sufficiently large. 
Hence to estimate |Ti — Ti| we divide the integration domain [0, t] into three parts, 
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[0, Ts ] , [Ta , t — Ts ] , and [t — , t] . As the estimates on different intervals are very 
similar, we only consider the first interval. 

h :=Z f ds[K(t-s)-K(s)]Re(W,e 1 ^W) f ds 1 q Sl [1 - X t(si)} 

JO J s 

pT'S _ _ pt 

-Z / - s) - #(s)]M(s) / d ai - xt(si)} . 

JO Js 

By (3.11) we have that 

\K(t -s)- K(t)\ <H(t - sy 1 *^^ — r- + (t - s)- 1 - r 1 + H 

<ri(l+«), 

because s < Ts, and t > T. This, together with the fact that |if (t — s) — ^(£)| < 
t~2g, implies 

\K(t-s) - K(t)\\Re(W ,e^W}\ + \K(t - s) - K(t)\\M(s)\ < Hs - a . 

Plugging this into Ji we obtain that 

i 

/■T3 

\h\Zt~ 1 - 6 / S -2|| g .|| (SiT = t-i-'52T5||g.|| (5iT <t-|- <5 T-5||g.|| 5iT . 
All together, we conclude that 

where e(T) — > 0, as T — >■ oo. 

This completes the outline of our proof. □ 
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